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SUMMARY 


Experimentally determined mass- size distributions of fragments ejected 
from craters formed in rock by hypervelocity impact have been combined with 
estimates of the rate of impact and mass distribution of interplanetary debris 
which strikes the lunar surface to find the rate and mass of fragments sprayed 
up from the lunar surface. It is shown that the flux of particles of a given 
mass ejected from the lunar surface will be at least three and probably four 
orders of magnitude greater than the flux of the impacting interplanetary debris 
of the same mass. Experimentally determined distributions of mass with velocity 
indicate that almost all of the debris is ejected at less than lunar escape 
velocity (2.4 km/sec) and contributes to secondary impact events. A small 
fraction of the ejecta, however, will leave the lunar surface at velocities 
greater than the escape speed. 

These results imply the presence of a lunar dust cloud of flying 
particles. The major fraction of the cloud is estimated to be a few kilo- 
meters deep with a spatial density at the lunar surface of the order of 10 s to 
10 7 ” times the spatial density of the interplanetary debris. These particles, 
together with the interplanetary debris, provide a powerful eroding agent that 
will continually abrade the lunar surface and reduce rubble to finer sizes. The 
inference to be drawn is that there must be at least a thin layer of dust-sized 
particles on the lunar surface — a layer which is being constantly agitated and 
stirred, by the impacts to form a heterogeneous mixture of material from the 
entire surface of the moon. It is to be expected that lunar probes and vehicles 
will be subjected to higher rates of impact when immersed in the dust cloud than 
during the earth-moon voyage . 


INTRODUCTION 


The moon, devoid of any appreciable atmosphere (ref. l) , is naked before a 
continuous bombardment by interplanetary debris of all sizes. The craters on 
the lunar surface visible to the limits of telescopic resolution are mute 
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testimony of this bombardment by large bodies over geologic time. Opik (ref. 2) 
and, most recently. Shoemaker and Hackman (ref. 3 ) have demonstrated that the 
number of what are believed to be post-mare impact craters increases exponen- 
tially with decreasing size, and. Shoemaker and Hackman have shown that the size- 
frequency distribution is in fair agreement with the distribution for terrestrial 
craters of apparent meteoric origin. It is significant that the correlation 
between lunar and terrestrial craters depends on the mass-frequency distribution 
of asteroidal bodies and recovered meteorite falls (refs. 4 , 5 ), for on this 
basis, by extrapolation to smaller bodies and considering the higher flux rates 
observed for meteors and micrometeoroids (e.g., refs. 6 , 7 and 8 ), it is to be 
expected that the craters and pits formed by smaller bodies and microparticles 
are the primary cause of the "microscale" relief and texture of the lunar sur- 
face (refs. 9 , 10). Indeed, the high flux rates for the smallest particles 
unimpeded by an atmosphere, probably constitute one of the principal, if not the 
primary, agents of erosion on the lunar surface. Meteoroid impact is, therefore, 
one of the major factors to be considered in the evolution of the moon to its 
present state, and a thorough understanding of the phenomena related to impact 
is essential to the study of the lunar surface. 

As part of a general study of the mechanics of hypervelocity impact, Ames 
Research Center has undertaken an investigation of meteoroid impact in a coop- 
erative program of research with the U. S. Geological Survey. Particular 
emphasis is being placed, on applications to lunar problems, specifically as 
applied to the NASA's programs for manned and unmanned exploration of the moon. 
Projectiles, accelerated to velocities up to 7 km/ sec in a nitrogen atmosphere 
of 50 to 75 mm of Hg, have been fired into various types of rock and sand targets 
using the light-gas gun facilities of the Ames Hyper velocity Ballistic Range. 

Some early results of this investigation have been reported by Shoemaker, et al. 
(ref. ll), and, Moore, et al. (refs. 12, 13, and 14). 

An important phase of the experimental program has been obtaining high- 
speed camera records of the ejection of material, from an impact crater. Even 
a cursory examination of the photographic records impresses the viewer with the 
amount of debris ejected from a crater formed in rock. From the instant the 
projectile contacts the target block, a cloud of rock flour and progressively 
larger fragments are spewed out of the embryonic crater until the crater reaches 
its final dimensions. The photographic records vividly demonstrate that for each 
lunar impact event a tremendous number of fragments should be ejected into bal- 
listic trajectories across the lunar surface. With a sufficiently high flux rate 
of interplanetary debris, one can envision a mass of fragments above the lunar 
surface which, continuously replenished, by ejecta from subsequent impact events, 
would approach a steady-state cloud of impact debris. The spatial density of the 
ejected particles should be expected to be considerably greater than that for 
the interplanetary debris, particularly near the lunar surface, and. the secondary 
impacts by the fragments should be many times more numerous than those by the 
interplanetary bodies. 

It is the purpose of this report to present an analysis of the spray 
fragments ejected, from the lunar surface by meteoroid impact based, on data avail- 
able from the joint NASA-USGS impact studies. Ejecta fragments collected from 
impact in rocks have been analyzed to ascertain the mass- size distribution of the 
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fragments. This analysis has been combined with estimates of the mass 
distributions of the solid material which impacts the lunar surface to give an 
order-of -magnitude estimate of the flux of the ejected fragments. Further, a 
study of ejected mass using the high-speed camera records has produced estimates 
for the mass-velocity distribution and velocity-ejection angle distribution of 
the fragments and, hence, of the spatial density and mass distribution of the 
ejecta across the lunar surface . 

The symbols used in the text are defined in appendix A. 


FLUX OF THE EJECTED FRAGMENTS 
Analytical Relationships 


It will be convenient first to express in analytical form the flux of 
fragments ejected from impact craters formed by a steady bombardment of projec- 
tiles. Empirical relationships will then be introduced to evaluate the flux of 
fragments from the lunar surface. Toward this end the functional forms 

m = M e g 1 (m e ) (l) 

Me = g 2 <>p) (2) 

dp ~ Ss( m p) (3) 

are introduced where mg ' is the mass of an ejected fragment, m is the 
cumulative mass of the fragments with individual masses equal to or less than 
mg. Mg is the total mass ejected from a crater, mp is the mass of a projec- 
tile, and fp is the cumulative number of projectiles (flux per unit area and 
time) with masses equal to or greater than Dip. Equation (l) implies that the 
fragmented products of an impact may be described by a simple comminution law 
independent of the crater dimensions. The use of this approximation for the 
present analysis will be discussed more fully when the empirical form of 
equation (l) is introduced. The second equation is based on empirical data 
obtained from impact-, chemical-, and nuclear-explosion cratering experiments 
which suggest Eg can be related to expended energy. If it is assumed that 
there is a characteristic impact velocity (constant specific kinetic energy) then 
Mg is dependent only on the projectile mass. The last, equation (3), is the 
form employed in the past by most investigators to describe the meteor or mete- 
oroid population in terms of particle mass, nip (e.g., refs. 6, 7, 8, and 15). 

Now if N is defined as the cumulative number of fragments with masses 
equal to or greater than 1% from any one crater, then the number of fragments 
of mass me from any one crater becomes -dN/dmg. The number of fragments of 
mass me from any crater can also be derived from equation (l) to be 
(l/me) (dm/dme ) . By equating the two expressions, integrating, and making use 
of equation (2) 
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N = g 2 (m p ) 


,m b gl « (me) 


me 




dffle 




= the cumulative number of fragments of ejecta 
with masses equal to or greater than 
formed as the result of one impact by a 
projectile of mass mtp. 

The quantify represents the mass of the largest fragment ejected from the 

crater, and the prime denotes differentiation with respect to the independent 
variable mg • 

The n umb er of impacts P per unit area and time contributed by . projectiles 
of a given mass m p is obtained by differentiating the cumulative number of the 
impacting projectiles, equation (3)- Thus 

df-p . . 

p = - s; = -ss’h?) 

and. one may write 

PN = -g 3 * (m p )g 2 (mp) f -Sl - i - Se) - dmg ( 5 ) 


= the cumulative number of ejecta fragments per 

unit area and time with masses equal to or greater 
than nig from all craters resulting from the 
impact of projectiles of mass nip. 

Finally, when f e is defined, as the cumulative number of all ejecta fragments 
per unit area and. time with masses equal to or greater than mp, 



and. 


fe = - 


g3’( m p)Ss( m p) 


,mb gi'(me) 


V 


mg 


mg 


dmedm-p 


( 6 ) 
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Equation (6) is a solution for the flux of fragments with masses equal to or 
greater than which are ejected from craters formed by a steady flux of 

projectiles with masses equal to or greater than mtp . It is to be noted that 
mb will depend on the size of the crater; that is, 

( 7 ) 


“b = S 4 (Me) 


and by equation ( 2 ) the upper limit for the first integration can be expressed 
in the terms of m p . In addition, the lower limit for the second integral 
requires special consideration. No values of Bn, which produce maximum fragments 
b^ SBialler than mg, may be included in the integration. Therefore, when 
equations ( 2 ) and (7) are coBibined to express m^ as a function of an, 
must be set equal to and the resulting expression for nip used as the 

lower limit of integration in the outer integral of equation (6). 


Empirical Relationships 


Material ejected from craters formed in rocks by the impact of hypervelocity 
projectiles has been collected and sieved through standard screens (e.g., ref. 1 6 ) 
to determine the Biass-size 2 distribution of the fragments. Figure 1 is a photo- 
graph of the ejecta from a typical crater formed in basalt with the fragments 
separated into logarithmic sized intervals ranging from large spall plates to 
finely ground rock flour. Experimental mass-size distributions of ejecta from 
eight craters are presented in logarithmic form in figure 2 as histograms of the 
fractional contribution from each size class and in terms (circle symbols) of the 
cumulative mass of fragments smaller than a given fragment size e. For both 
presentations departures from linearity appear to be caused by erratic spalling 
of the larger fragments and by preferential fracturing along grain boundaries for 
the smaller fragments (ref. 12). Analysis of the debris finer than approximately 
50 microns has not been completed, but results obtained to date suggest that 
there is a cutoff in particle size near 0.1 micron and a consequent departure 
from the logarithmic linearity in the sizes less than 30 microns (fig. 3). The 
apparent cutoff at 0.1 micron, however, may be attributable partly to incomplete 
collection of fine material. 


Results similar to those given in figures 2 and 3 have been reported by 
Roberts and Blaylock (ref. 17) for the ejecta recovered from the Stagecoach series 
of high explosive cratering experiments in desert alluvium. Moreover, some 
unpublished results obtained by Shoemaker indicate that, for explosive debris in 
a pumiceous tuff, a linear logarithmic presentation is valid well into the Biicron 
region. To a good approximation, adequate for this analysis, the mass-size 
distribution of cratering events over a considerable range of expended energies 
Biay be expressed (normalized with respect to the size b of the maximum fragment) 

/ \<x 


m = Me 



( 8 ) 


2 The size of a fragment is taken to be the intermediate dimension of the 
particle. 
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■with m and. Mg, as before, the cumulative and total masses ejected, e the 
fragment size, and a a constant. Equation (8) is the usual form employed to 
describe approximately the comminution of brittle, macroscopically homogeneous 
solids in commercial or laboratory crushing mills (e.g., ref. l8). When similar 
geometry for all fragments is assumed, equation (8) can be rewritten 



with e and b being replaced by their respective masses, mg and m^ as required 
for solution of equation (6) . 

A relationship between total ejected mass M e and expended energy is shown 
in figure 4 which presents results from both laboratory impacts and shallow 
subsurface explosive cratering experiments in various media (unpublished data 
by Moore and refs. 19-30) . The excavated mass for explosives is relatively 
insensitive to the physical properties of the medium but the charge burial depth 
is an important factor. There is no direct analogy between buried explosive and 
impact cratering events, but the effective scaled depth of burial A for 
explosives to agree with impact data appears to be quite small (A = 0 to 0.25) 
based on the total mass ejected. 

Theoretical considerations (e.g., ref. 3l) and experimental evidence (e.g., 
refs. 32 and 33) suggest that for impacts the excavated mass Me should be 
directly proportional to the expended energy as shown by figure 4 . However, 
because of the high specific energies that occur during impact at meteoric veloc- 
ities, comparable to and perhaps exceeding those for the nuclear experiments, 
considerably more energy will be irreversibly trapped as heat and lost to the 
cratering processes than occurred in the events presented in figure 4 . Moreover, 
a major fraction of the available energy from an impact is lost as kinetic energy 
in the ejecta squirted out at velocities up to three times the projectile veloc- 
ity. For these reasons, some deviation from a direct proportionality between 
energy and ejected mass might be expected for impacts at meteoric velocities so 
the form adopted for equation (2) is 


Mg = kfmp Q Vi 2 ) = g 2 (mp) (10) 

with k-L and (3 being constants and Vi the impact velocity. 

The expression for g (mp) , as indicated previously, will be taken in the 
form commonly used to describe the flux of meteors and meteoroids (refs. 6, 7 , 8, 
and 15) • 

fp - kgHipA = g 3 (m p ) (ll) 
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The final functional relationship g 4 (Me) can he Expressed 

mb = ksMg = g 4 (Me) ( 12 ) 

based on the results of a correlation presented as figure 5 (unpublished data 
by one of the authors (Moore) and refs. 27-30, 3^, and 35 ). In equations (ll) 
and ( 12 ) k 2 , k 3 , 7 , and & are constants. 


Numerical Results 


Inserting equations ( 9 ) through (12) in equation ( 6 ) and carrying out the 
indicated integrations, one obtains for the flux of impact- spray fra g ments of 
mass equal to or greater than nig (see appendix B) 

f e = k 2 Kme £+7 (13) 


with 




1 

.7 + 1 - ( 08 / 3 ) 


1 

7 + 1-5, 


€ 


7+1 

5 


1-7 


for 

7 + 1 - — < 0, 7 + 1 - S < 0, and a / 3 

Since fp = k 2 mp 7 , equation ( 13 ) can be recast in terms of the frequency of 
the impacting projectiles 


f e = K 


me 


e+7 




m, 


L P 


and when projectiles and ejecta fragments of the same mass, mp e , are considered, 
^e/fp ~ K(mp e ) 

cumulative number of fragments with masses > 

= (lk) 

cumulative number of projectiles with masses > mp e 

for unit area and time and for 


mg = m p = m pe 
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Seven constants (a, (3, 7 , S, k x , and Vi) and the independent variable 
ittpe are required to evaluate equation (l4) . Values for a range from 0.3 to 
0.7 for the laboratory craters in basalt (figs. 2 and 3) and for the explosive 
craters in alluvium (ref. 17)* The ratio f e /fp is relatively insensitive to 
d in this range and a value of a of 0.4 will be adopted. The quantity 3, for 
reasons already cited, should be less than 1.0 and for large events may be lower 
than 0.88 which is suggested, by empirical scaling laws for explosion craters in 
alluvium (ref. 36). For the present calculations, however, a value of 1.0 will 
be used in conjunction with a conservative value for k 1 (see fig. 4) 

k x = — - — gm/erg 

8xi0 8 

The literature on meteor and meteoroid frequencies indicates that 7, k 2 , 
and Y± will depend on the type of projectiles being considered. The geocentric 
flux suggested by Hawkins (ref. 37) and Brown (ref. 4) for debris of suspected 
asteroidal origin is shown in figure 6 together with those of Whipple (ref. 6), 
Hawk in s and Upton (ref. 7)> and McCracken and Alexander (ref. 8) for the smaller 
particles usually associated with cometary and interplanetary dust. Recent 
results reported by Soberman and Hemenway (ref. 15) which pertain to microparticle 
collection with the Venus flytrap are also included for comparison assuming spher- 
ical particles with a density of 3 gm/cm 3 . The values of 7 and k 2 for the fp 
distributions shown in figure 6 are tabulated below 


Type of debris 

Author 

7 

l0 Sl0 k 2 

Meteorites 

Brown 

-0.8 

-16.1 

Asteroids and fireballs 

Hawkins 

-1.0 

-14.5 

Meteors 

Hawkins and Upton 

-1.34 

-12.2 

Meteors 

Whipple 

-1.0 

-11.3 

Interplanetary dust 

McCracken and 

-1.7 

-17.0 


Alexander 


-1.6 

Micrometeoroids 

Soberman and 

-0.4 


Hemenway 

These values represent geocentric or near earth populations and the question 
arises concerning the applicability to the lunar environment. Whipple (ref. 38) 
has concluded that there is a geocentric concentration of the smaller particles 
and that the number of particles decreases as the 1.4 power of the distance from 
earth. At the distance of the moon from earth, Whipple finds the population of 
particles with masses larger than 10“ 9 grams would be about the same as that esti- 
mated for the zodiacal dust cloud, approximately three orders of magnitude less 
than the near earth values. Whipple's interpretation, however, has been questioned 
by McCracken and Alexander (ref. 8) on the basis that there is no reason to sus- 
pect a priori that the mass -magnitude (visual or photographic) relationship is 
constant (7 = -1.0) over a large range of particle masses. The measurements of 
particle frequencies reported by McCracken and Alexander from satellites, probes, 
and sounding rockets yield 7 = -1.7 directly. Hawkins and Upton's analysis also 
suggests that 7 departs from -1.0 in support of the argument of McCracken and 
Alexander . 
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For the present analysis, fortunately, the absolute populations are not 
required (i.e., the value of k 2 ) to evaluate the ratio fe/fp, but only the mass 
distributions (i.e., /) are required. The view adopted herein is that, regardless 
of any geocentric concentration, the mass distribution in the lunar environment 
should be essentially unchanged from near earth distributions, at least to parti- 
cles of 10" 10 grams. For smaller masses, the terrestrial atmosphere may contrib- 
ute to a difference between the terrestrial and lunar distributions. The results 
presented by Soberman and Hemenway would seem to fall in this category and, hence, 
will not be employed herein. However, since the ratio fe/fp is quite sensitive 
to 7 and Vi, the ratio will be evaluated for the remaining values for 7 in 
combination with the appropriate modal velocities Vi for asteroidal and cometary 
debris. 

Altho ugh population distributions are assumed unchanged, consideration must 
be given to the difference in gravitational attraction between the earth and moon 
as it affects the impact velocity Vi- Geocentric velocities for asteroidal debris 
may be taken to be 1 9 km/sec. While the velocities for cometary debris may be 
as high as 72 km/sec, a value of 30 km/sec will be employed herein as represent- 
ative for most cometary material. 

From the conservation of energy 

Vj_ = 4 2gr + V TO 2 

where g is the gravitational acceleration at the surface of a body with a 
radius r, and is the approach velocity at an infinite distance from the 

body. The geocentric velocities, therefore, correspond to approach velocities 
of approximately 15 km/sec for the asteroidal debris and 28 km/sec for represent- 
ative cometary material. The gravitational field of the moon contributes so 
small a velocity increment that impact velocities of 15 and 28 km/sec are ade- 
quate for use in equation (l4). 

One final set of constants 8 and k3 is required. The correlation presented 
as figure 5 indicates that 8 has a value varying from about 1.0 for the smallest 
craters to 0.8 for the largest craters. A value of 1.0 will be employed for the 
numerical evaluations of equation (l4) together with k3 = 10" 1 for the cometary 
debris. Values of 5 = 0.8 and k 3 = 2X10" 1 will be used for the asteroidal com- 
ponents . It is to be noted that a change in the relation of to Mg should be 
expected for the larger craters (Mg > 10 4 grams; m.p > 10 grams) because of the 
probability of secondary breakup along pre-existing fractures or joints in the 
country rock, and the change of the importance of material strength with the size 
of ejected fragments. 


A summary of the numerical reduction of equation 
following table : 

( 14 ) is 
Vi, 

listed in the 


Debris 

Source 

7 

km/ sec 

K 

e 

Asteroids and fireballs 

Hawkins 

-1.0 

15 

1.8x10 3 

0 

Meteors 

Whipple 

-1.0 

28 

4.9x10 s 

0 

Meteors 

Hawkins and Upton 

- 1.34 

28 

1 .JXL 0 4 

0 

Interplanetary dust 

McCracken and 
Alexander 

- 1.7 

28 

6.oxio 4 

0 
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and the results are presented In graphical form as figure 7* Wo values are 
shown for Brown's distribution (7 = -0.8) because the inequality (appendix B) 

7 + 1 - (ab/3) <0 is violated with the constants required to fit this case . 

The flux of the ejecta fragments is, as expected, greater in all cases than 
the flux of the impacting interplanetary debris . Because of the choice of 
constants, the exponent e has a value of zero and the ratio f e /fp for -each 
distribution is independent of the mass nip e . 

Hawkins’ distribution for asteroidal projectiles gives a ratio of f e /fp 
somewhat greater than 10 s . Although the extrapolation of the small-scale labora- 
tory experiments to the scale of the telescopically visible lunar craters is 
tenuous, the numerical results clearly suggest that the lunar surface is littered 
with large fragments thrown out of the impact craters. The great bulk of the 
larger fragments will, as will become evident, lie on the flanks of the craters 
but some can be expected at considerable distances from the points of impact. 
Secondary breakup, as mentioned previously, should be an important factor for 
large craters and it is to be expected that the ratio f e /fp should tend toward 
zero for the largest lunar craters (increasing mp e ) . 

For the smaller projectiles associated with cometary debris, the distribution 
of McCracken and Alexander gives an ejecta flux f e of almost five orders of 
magnitude greater than the cumulative projectile flux fp. The ratio f e /fp for 
the Hawkins and Upton distribution is approximately 10 4 while, as a lower limit 
for the microparticles, Whipple's distribution gives f e larger than fp *y 
between three to four orders of magnitude. Thus, in all cases the secondary 
fragments exceed the projectile population by between 10 3 and. 10 4 . Since most 
of these fragments would be dustlike particles, a lunar dust cloud might be pro- 
duced by a sufficiently high flux rate of projectiles. 


DISTRIBUTION OF EJECTA FRAGMENTS 
ACROSS THE LUNAR SURFACE 

Trajectories of Ejecta Fragments 


In the absence of all but lunar gravitational forces, a fragment that is 
ejected from the moon and returns to the lunar surface will travel along a 
segment of an elliptic orbit. Except for certain limiting cases the perilune of 
the orbit will be within the body of the moon, and apolune will correspond to the 
position for maximum altitude of the fragment above the lunar surface . The 
circumferential distance or range R from ejection to the secondary impact is 
(see, e.g., ref. 39) 


R = 2r tan -1 



sin 8 cos 9 


V e 2 cos 2 0 


(15) 
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where 


0 <! 0 <? jc/2 


and 



Ve 2 

rg 0 


V e is the ejection velocity, g Q the gravitational acceleration at the lunar 
surface, r the lunar radius, and 0 the angle of ejection measured with respect 
to the local horizontal. The maximum altitude h^x may be expressed 


h. 


■max 


_ r 


v e 2 - 1 + 

”11/2' 

1 - V e 2 (2 - V e 2 ) cos 2 0l 

■ 

2 - V e 2 


(16) 


Since r and g Q are known (l.74xi0 3 km and 1.62X10 -3 km/sec 2 ), the determination 
of the range and maximum altitude of a fragment depends only on knowledge of the 
vector velocity of the ejection from a crater. It is to be noted that, indepen- 
dent of 6 , fragments are lost to the gravitational field when 

V e = \/2rg 0 = 2.38 km/sec 

As mentioned in the Introduction, high-speed camera records have been 
obtained of the rock fragments spewed out of impact craters and have provided a 
basis for estimating the mass-velocity distribution of the ejecta. The photo- 
graphic records have been obtained at three nominal framing rates (lO 4 , 10 5 , 
and 10 s frames/sec) and a sequence at the highest rate is shown as figure 8. 
Although the total time increment covered by the record is approximately 70 micro 
seconds, the final crater geometry in the basalt target probably has been estab- 
lished and only the dispersion of the larger fragments moving at relatively low 
speeds remains to complete the ejection sequence. 

Results from an analysis of the photographic records are presented in 
figures 9 and 10 for aluminum spheres (l/l6- and l/8-inch diameter) impacting at 
normal incidence on basalt at a nominal velocity of V± = 6.25 km/sec. Figure 9 
presents the ejection velocity V e as a function of the ejection angle 0, and 
figure 10 presents the cumulative mass m which is ejected with velocities in 
excess of V e . Material ejected at the highest velocities (V e ~ 19 km/sec) and 
lowest angles (0 < 20°) is believed to be associated with a jet emanating from 
the interface between the projectile and target (see, e.g., refs. 40, 4l, and 42) 
Subsequent to the jetting the ejection velocity quickly decays to very small 
values and the main mass of fragments leaves the craters with velocities less 
than 0.5 km/sec and with 0 greater than 45°, The discontinuous variation in 
the cumulative mass with velocity is believed to be associated with a transition 
from plastic to elastic flow behind the shock front which propagates outward from 
the point of impact. 
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Estimates for the distribution of mass with velocity for the 15 and 28 
km/sec impact velocities are also included in figures 9 and 10 . In extrapolating 
the laboratory results to higher impact velocities, it is convenient to assume 
that the expended energy remained constant (l/ 2 mpVi 2 ~ 9 >< 10 s erg). Thus the total 
ejected mass Me becomes independent of the impact velocity (eg,. (10)) and the 
projectile mass m.p can be considered to vary inversely with the square of the 
impact velocity Vi. Moreover, since the main bulk of the debris is ejected 
under conditions of elastic flow, the absolute values for the ejection velocity 
V e were considered essentially independent of the impact velocity except for a 
small fraction of the ejected mass associated with jetting and plastic flow. It 
is to be expected that, since jetting is a hydrodynamic process, the velocity and 
mass of the jetted material will scale approximately linearly with some charac- 
teristic velocity and the cube of a characteristic length (i.e., projectile mass), 
respectively. Thus, the maximum ejection velocities for 15 and 28 km/sec in 
figures 9 and 10 have been scaled by taking V e max = (Vi/6.25) while the associ- 
ated jetted masses in figure 10 have been scaled by m m -in = (6.25/Vi) 2 . The cal- 
culated end points for jetting and the experimentally determined curves have 
been joined arbitrarily to complete the velocity-ejection angle relationships and 
the cumulative mass-velocity distributions for the higher impact velocities. 

It should be cautioned that the above procedure neglects any effects of the 
increase in specific kinetic energy for the higher impact velocities. An increase 
in impact velocity increases the fraction of the kinetic energy which is trapped 
irreversibly as heat during the shock compression of the projectile and target 
media. Thus, even though the kinetic energy may remain constant, the fraction of 
energy available for crushing and ejection of fragmented material decreases with 
an increase in impact velocity and one should expect a slight decrease in the 
total ejected mass M e . The effect would be to translate the estimated mass- 
velocity distributions shown in figure 10 to the left and reduce the cumulative 
masses which are ejected at velocities in excess of any given ejection velocity. 
The over-all change is believed to be small, but it should be kept in mind that 
results presented herein for 15 and 28 km/sec may represent the maximum mass 
which can be ejected at any given velocity. 


Distribution of Ejecta Fragments 
in the Lunar Environment 


The circumferential distance ejected fragments will travel across the moon 
from an impact crater (derived from the estimated variation of cumulative mass 
with vector velocity and from eq. (15)), is presented in figure 11 as the normal- 
ized cumulative mass m/mp that contributes to secondary impact events beyond 
a given distance from the point of projectile impact. Similarly, with the use of 
the experimental data and equation (16) , the apolune distance of the ejected frag- 
ments is presented in figure 12 as the normalized cumulative mass m/mp of the 
ejected fragments ejected above a given lunar altitude. Three curves are shown 
in the presentations, one derived directly from the experimental data for a nomi- 
nal impact velocity of 6.4 km/sec and estimates for 15 km/sec and 28 km/sec 
impact velocities. 
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The results presented in figures 11 and 12 indicate that the main fraction 
of the ejecta, say 90 percent, will never exceed an altitude of the order of 10 
kilometers and will return to the lunar surface within a radius of the order of 
30 kilometers from the point of impact. If one considers virtually all of the 
ejecta mass (99 percent) , the maximum altitude and radial distance increase to 
30 and 150 kilometers, respectively. 

The flux of the ejected fragments has been found to he from 10 s to almost 
10 5 times the flux of the interplanetary debris (fig. 7 ) . Since the velocities 
of the bulk of the ejected fragments are some two orders of magnitude less than 
the velocities of the projectiles, it follows that the spatial density of the 
ejecta at the lunar surface is from 10 5 to 10 7 times the spatial density of the 
impacting interplanetary debris. 

The last 1 percent of the ejected mass includes perhaps the most significant 
fraction of the fragments. Dependent on the impact velocity, a mass of fragments 
equivalent to several projectile masses is ejected at velocities in excess of 
lunar escape velocity and is lost to the lunar gravitational field. It is inter- 
esting to note in this connection that the impact experiments and analysis herein 
for impact into solid basalt indicate that the velocity for which the mass lost 
from the moon is equal to the acquired projectile mass is between 9 and 10 km/sec. 
Thus, it appears that all asteroidal and cometary projectiles are effective in 
removing mass from the moon. 


DISCUSSION 


The results from the present study indicate the existence of a lunar 
"atmosphere" of high-speed ejecta. Moreover, the results suggest that the lunar 
surface material is composed of a mixed rubble of unsorted rock fragments ranging 
in size from large blocks to submicron sized particles. A model of the lunar 
surface in which a gravelly material extends to depths of at least several meters 
has been proposed recently by Baldwin (ref. 43) based on observations of the 
thermal radiation from the moon. In view of the higher flux rates for the 
smallest projectiles, both interplanetary solids and crater ejecta, it is to be 
expected that the surface rubble is continually abraded to finer sizes and that 
over geologic time small voids and cracks have gradually filled with pulverized 
rock. Indeed the abrasive action of the "atmosphere” must contribute, again 
consistent with optical and radio telescope observations, to at least a thin 
mantle of fine particles over the entire lunar surface--a layer which is con- 
stantly agitated and stirred by impact to form a heterogeneous mixture of mare 
and terra material. 

A question immediately arises concerning the applicability of the current 
analysis, based on impact in solid targets, to impact in an unconsolidated or 
weakly cemented particulate medium. Some results for craters formed in a weakly 
bonded quartz sand (crushing strength approximately 10 7 dynes/cm 2 ) are shown in 
figure 4 by the filled symbols. Sand particles with a grain size of 0.3 mm were 
employed for these tests in combination with projectiles which had a diameter of 
l/8 inch (3.18 mm). For a given expended energy, the mass Me ejected from a 
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crater in the bonded, sand exceeds the mass from the craters in basalt by a factor 
of at least three. In addition, the material sprayed out from the craters in 
sand contains few, if any, large spalls and consists almost entirely of finely 
crushed quartz, disaggregated sand grains, and a few clusters of sand grains 
welded together with fused material (silica). Both the increased mass Me and 
the lack of spalls will serve to increase the calculated values of f e /fp over 
those based on the solid targets. The increase would be about an order of magni- 
tude for McCracken and Alexander’s distribution (7 = -1.7) and a factor of, per- 
haps, four for Whipple 1 s distribution (/ = -1.0). On this basis the numerical 
values for fe/fp in figure 7 should be conservative. 

There is, however, no assurance that the data presented in figure 4 for 
impact in sand are. valid when the projectile is small relative to the grain size 
of the target and/or the target material is a diffuse low-density medium, such as 
that suggested by D. D. Cudaback (ref. 44) for the lunar surface. It is probable 
that for these latter conditions less mass and fewer fragments will be ejected by 
an impact as compared to an equivalent event in a solid rock. This would tend, 
of course, to reduce the effectiveness of the smaller particles of interplanetary 
debris in producing secondary fragments. Any reduction should be offset by the 
increases noted in the previous paragraph but a quantitative evaluation of the 
combined effects integrated over a spectrum of projectile masses cannot be made 
at the present time. It is worth noting, however, that the ratio of projectile 
to sand-grain dimensions for the data in figure 4 is the same as that which would 
occur with a lunar dust layer consisting of 1 micron particles and 10 9 gram 
interplanetary projectiles (10 micron diameter, 3 g/cc density). 

One difference observed between impact in rock and sand is a change in the 
velocity-ejection angle distribution of the fragments from that shown in figure 9. 
Based on the results currently available for impact in sand, the ejection angle 
d increases rapidly from the low values associated with jetting, attains a peak 
value, and then decreases monotonically to some minimum value. Such a variation 
in 9 will increase the range R and decrease the apolune height h^x from 
that shown in figures 11 and 12 for an undetermined fraction of the low- speed 
ejecta. This difference in the ejection angles, attributable to reduced target 
strength, does not introduce any important changes to the analysis and, in parti- 
cular, will not change the previous estimates of the ratio f e /fp or the spatial 
density of the ejecta fragments at the lunar surface. 

A lunar atmosphere of high-speed projectiles poses a hostile environment 
for both manned and unmanned exploration of the moon. The frequency of impacts 
on manned vehicles and instrumented probes will be far greater on the moon than 
that encountered during the earth-moon voyage. The secondary impacts should not 
greatly increase the hazard from catastrophic punctures, because the bulk of the 
ejected material travels at low speeds, but the abrasive action of the multiple 
impacts, even at low speed, may present operational and maintenance problems for 
optical and solar-cell surfaces. Multiple impacts will steadily degrade such 
surfaces, and with sufficiently high rates of impact, may appreciably shorten the 
useful lifetime of the related equipment and experiments. Unfortunately, an 
estimate of the degradation cannot be made with any confidence at the present time 
because of the uncertainties in the absolute values for the flux of interplanetary 
solids onto the moon. With reference to figure 7, this can be demonstrated by 



Whipple’s flux distribution (ref. 6) and interpretation of satellite 
micrometeoroid measurements (ref. 20 ) ; one obtains an absolute flux and a spatial 
density for secondary particles having masses greater than 10" 10 grams of the 
order of one particle per second over an area of approximately 10 square meters 
and one particle per 10 s cubic meters, respectively. In contrast, the direct 
application of McCracken and Alexander's results (ref. 8) yields an impact rate 
and spatial density almost six orders of magnitude higher — 10 impacts per second 
per square centimeter and one particle per 10 3 cubic centimeters, repectively. 

The lowest rate might be considered minor as compared to the highest rate. The 
actual values are probably somewhere between the two extreme s , 3 but it remains 
for impact counters in the lunar environment to establish the rate at which ejecta 
are thrown up by the impact of interplanetary solids. It is, perhaps, superfluous 
to note that the investigation of impact on and near the lunar surface is a prereq.- 
uisite to the undertaking of manned exploration of the moon. 


CONCLUDING REMARKS 


The analysis of the distribution of solid material thrown out of impact 
craters on the moon has ramifications which extend well beyond the scope of the 
present paper. Nevertheless, it should be mentioned in closing that the impli- 
cations arising from the indicated negative accretion of mass by the moon are 
particularly significant. An effectively steady rate of ejection of an appre- 
ciable mass of solids from the lunar environment over geologic time is fundamen- 
tally important both as an evolutionary process which has brought the moon to its 
present state and as an origin for some of the solid debris in interplanetary 
space. The latter, of course, has a direct bearing on many problems, notably the 
zodiacal dust, the gegenschein, the Kordelewski clouds, and any geocentric concen- 
tration of dust. Moreover, the possibility of ejecting several projectile masses 
of material away from the gravitational field of the moon provides an origin for 
some meteorites. It is interesting to note, also, that impact appears to be an 
extremely efficient mechanism for the levitation of small particles above the 
lunar surface which is necessary in Gold’s theory for a space-charge fluidization 
and transport of dust on the moon. These and other related problems serve to 
emphasize the role of lunar impact as one of the important keys to the exploration 
and better understanding of the solar system. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Dec. 27, 1962 


3 <John A. O’Keefe, NASA Goddard Space Flight Center, has pointed out (private 
communication) that the light scattered by a dust cloud corresponding to the 
highest impact rate should provide a readily detectable extension of the lunar 
cusps near first quarter. No such extension has been observed. Refined photo- 
metric techniques should be capable of detecting a lunar dust cloud, some three to 
five orders of magnitude less dense than that corresponding to the highest flux 
rate estimated above. 
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APPENDIX A 


SYMBOL NOTATION 


The principal symbols employed throughout the body of the text are 
su mm arized as follows: 


To 

d 

e 

E 

fe 


So 

hmax 

K,ki,k 2 ,k3 

m 

in 

“b 

% 

mp 

^pe 

Me 

r 

R 

V e 

Vi 


size of largest fragment ejected from a crater 
explosive burial depth, ft 
size of fragment ejected from a crater 
expended energy expressed in pounds of TNT 

flux of fragments per unit area and time which have masses equal to 
or greater than me 

flux per unit area and time of interplanetary debris which have 
masses equal to or greater than mp 

gravitational acceleration at the surface of the moon 

maximum altitude a fragment attains above the lunar surface (apolune) 

constants introduced in equations (10) through (13) 

cumulative mass of fragments with masses equal to or less than 1 % 

cumulative mass of fragments ejected from a crater with velocities 
equal to or greater than Y e 

mass of largest fragment from a crater 

mass of fragment from a crater 

mass of projectile or interplanetary debris 

mass of fragment and projectile when 3% = mp = mp e 

mass of material excavated from a crater 

radius of moon or planet 

circumferential range of fragment ejected from a crater 
ejection velocity 
impact velocity 
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V velocity of the interplanetary debris at an infinite distance from 

00 

the moon or planet 
V e V e / V*io 

a, (3,7, 5, e constants introduced in equations (lO) through ( 13 ) 

Q ejection angle of fragment with respect to the local horizontal 

A scaled depth of explosion, d/E 1 ^ 3 

The cgs system is employed throughout except where noted. 
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APPENDIX B 


EVALUATION OF EQUATION (7) BY USE OF EMPIRICAL RELATIONS HIP S 


The flux of fragments ejected from craters formed by a steady flux of 
projectiles with mass co > mp !> mp i is expressed in equation ( 7 ) as 


f e " - 


mp. 


gs , (mp)g 2 (mp) P dmedmp (Bl) 

J me ^ 


and the empirical relationships which will be used in evaluating this equation 
(eqs. ( 9 ) through (12)) are 


/ \a/3 

- (t) 


g 2 (m p ) = Vj 2 ) = Me 

g 3 (m p ) = k 2 m p 7 = f p 


g 4 (Me) = mb = k s Me - k. 


k 


i m p 




(B2) 

(B3) 

(B4) 

(B5) 


Using equations (B2) to (b 4) in (Bl ) , one obtains 


fe 


a 




r r r mb ^ c “ /s) 


m r 


m 


L Pi 


Me 3 % 


(a/3) 


dm e dmp 


After the first integration 


“ 7*1*2 (| f V 


a - 3 


m Pl 


mbWa) - 1 . ^(a/s) - 1 


Mb 


(a/ 3 ) 


dm p (b6) 


provided that a f 3* For the data in this report and in reference 17, a is 
less than one. Equation (B5) can be used to recast equation (b 6 ) into 

00 

L Pi 



C 2 m p 


7-5 


- m-r 


7 - 


U.U 

3 


■= - 1 


Me 


dm. 


l P 


(B7) 
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with 



C 2 - 





( cc/3 ) - 1 


Since no projectiles may he included in the final integration which would 
produce fragments mb < equation (B 5 ) is rewritten in terms of me and mp 


mb - me - k 3 





so that 



(B 8 ) 


= mass of the projectile which produces a 
maximum fragment with = mg 

In tegrate on of equation (B7) using the equality of equation (b 8 ) in the lower 
limit gives 

f e = k 2 Kme £+7 (B9) 


with 




1 

-7 + 1 - (a6/3) 


1 

7 + 1-5- 


e 


7+1 

5 


1 


7 


for 

7+l--y<0, 7+1- S<0 

The two inequalities are a mathematical restriction that arises as a consequence 
of the form of expression employed to describe the flux of interplanetary debris 
(eq. (B4) ) and the final integration to bodies of infinite mass (eq. (B7))» This 
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can be illustrated by differentiating equation (b 4) to obtain the number of 
interplanetary projectiles P having a mass nip 


P 


dfp 

dmp 


- 7k 2 mp 


7-1 


(BIO) 


The over-all mass of the bodies whose individual masses are mp is 

7 

ratpP = -7k 2 mp 

and the total or integrated mass of all the bodies with masses greater than nip 
becomes 

f°° 7 

Total mass = -7k 2 / mp dmp (Bll) 

J mp 

When 7 < -I, a finite mass is obtained and with 7 > -1, an infinite mass 
results . 

Physically, equations (BIO) and (Bll) indicate that for 7 < -1, the number 
of particles and their mass increase without limit as mp tends to zero, while 
for 7 > -1, the number of bodies decreases but the masses involved increase 
without limit as mp tends to infinity. Thus, the latter implies (through 
eq. (B3)) an infinite mass excavated for the lunar impacts and an infinite number 
of secondary fragments. The specified inequalities for equation (B9) are directly 
related to this limitation in the mathematical development . The limitation can 
be circumvented, of course, by specifying a finite upper limit in equation (B7), 
but the added complexity is not justified for the present analysis. 
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Figure 1.- Ejecta from a typical crater formed in basalt with fragments separated into logarithmic- 
sized intervals. 
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Figure 2.- Mass-size and cumulative mass-size distributions of fragmented rock from craters formed 

basalt by hypervelocity impact . 












Cumulative mass of fragments smaller than 
Total ejected mass 


Round Sp - 445 



Figure 3.- Mass-size and cumulative mass-size distrib 
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Figure 4.- Correlation of the total ejected mass with expended energy for explosive and impact 

cratering events. 




Figure 5 . - Correlation of the mass of the largest fragment with total ejected mass for explosive and 

impact cratering events. 
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Figure 7*- The ratio of f e /fp for different types of projectiles. 
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Figure 8.- High-speed framing camera records of an impact in basalt; framing rate = 1.04X10 6 frames 

per second; impact velocity = 6.4 kilometers per second. 
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Figure 8.- High-speed framing camera records of an impact in basalt; framing rate = 1.04xi0 6 frames 
per second; impact velocity =6.4 kilometers per second - Continued. 
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Figure 8.- High-speed framing camera records of an impact in basalt; framing rate 1.04x10 frames 
per second; impact velocity =6.4 kilometers per second - Concluded. 
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Figure 9*~ Variation of ejection velocity with ejection angle for fragments 
thrown out of craters formed in basalt by hypervelocity impact. 





Figure 10.- Cumulative mass ejected in excess of a given velocity for fragments thrown out of craters 

formed in basalt by hypervelocity impact. 




Cumulative ejected mass 
projectile mass 



logio Range in kilometers, logio R 


Figure 11.- Distribution across the lunar surface of the mass ejected from impact craters formed 

in basalt. 
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Figure 12.- Distribution above the lunar surface of the mass ejected from impact craters formed in 

basalt . 



